Phase Separation in One-Dimensional Driven Diffusive Systems 
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A driven dilTusive model of three types of particles that exhibits phase separation on a ring is 
introduced. The dynamics is local and comprises nearest neighbor exchanges that conserve each of 
the three species. For the case in which the three densities are equal, it is shown that the model 
obeys detailed balance. The Hamiltonian governing the steady state distribution in this case is given 
and is found to have long range asymmetric interactions. The partition sum and bounds on some 
correlation functions are calculated analytically demonstrating phase separation. 
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Driven diffusive systems have been a subject of exten- 
sive studies in recent years Being driven by an 
external field, these systems are governed by dynamics 
which generically does not obey detailed balance, leading 
to steady states with non- vanishing currents. This class 
of systems provides a relatively simple framework within 
which collective phenomena far from thermal equilibrium 
may be studied |^,^. 

Theoretical studies of models of driven diffusive sys- 
tems reveal some basic differences between phase transi- 
tions taking place under equilibrium and non-equilibrium 
conditions. For example it is well known that phase tran- 
sitions and spontaneous symmetry breaking are not ex- 
pected to take place in one dimensional (Id) systems in 
thermal equilibrium at finite temperatures, as long as the 
interaction in the system is short range. Recently it has 
been demonstrated that this is not the case in systems far 
from thermal equilibrium j5|. Indeed, a simple example 
of an open system with non-conserved order parameter at 
the boundaries was shown to exhibit spontaneous sym- 
metry breaking [Q . 

A closely related problem is that of phase separation. 
An interesting question is whether Id homogeneous sys- 
tems (i.e. with no boundary effects as in a ring geometry) 
are capable of exhibiting phase separation, in cases where 
conserving dynamics is involved. It has been shown that 
inhomogeneities, such as defect sites or particles may 
trigger the formation of macroscopic regions of high den- 
sity bound to the defect |^ . Recent numerical studies of a 
model associated with sedimentation of colloidal crystals 
have suggested that phase separation may even occur in 
homogeneous non-equilibrium systems |Q . 

Phase separation is accompanied by coarsening phe- 
nomena in which the typical domain size grows indefi- 
nitely with time Examples of equilibrium Id systems 
with local dynamics which exhibit coarsening are the zero 
temperature limit of kinetic Ising models and the noise- 
less Landau-Ginzburg equation. However these systems 
do not coarsen at finite temperature. 

In the present Letter we introduce a simple model of 



phase separation in a Id driven diffusive system. This is 
a model of 3 species of particles on a ring in which near- 
est neighbors on the lattice are exchanged with specific 
rates. Thus the dynamics is local, fully stochastic and 
conserves each of the three species. It is shown that for 
the special case in which the average densities of the three 
species are equal, the dynamics obeys detailed balance. 
In this case the steady state distribution is shown to be 
given by a Hamiltonian which has long range asymmetric 
interactions. The phase separation which takes place in 
the model is explicitly demonstrated for this case. It is 
argued that phase separation takes place in the general 
case where the densities of the three species are unequal 
and that the typical domain size coarsens as Int. The 
model is easily generalized to n species and phase sepa- 
ration is found provided n > 2. 

A. Definition of the model: The model is defined on a 
Id lattice of length N with periodic boundary conditions. 
Each site is occupied by either an A, B, or C particle. 
The evolution is governed by random sequential dynam- 
ics defined as follows: at each time step two neighboring 
sites are chosen randomly and the particles of these sites 
are exchanged according to the following rates 



AB 
BC 
CA 



BA 
CB 
AC. 



(1) 



The rates are cyclic in A, B and C and conserve the 
number of particles of each type. 

Consider a system with Na particles of type A, Nb of 
type B, and Nc of type C. For q = 1 the particles un- 
dergo symmetric diffusion and the system is disordered. 
However for q ^ 1 the particle exchange rates are bi- 
ased. Since the model is invariant, for example, under 
the exchange A ^ B and q 1/q, it is sufficient to 
consider q < I. In this case the bias drives, say, an A 
particle to move to the left inside a B domain, and to 
the right inside a C domain. Therefore, starting with 
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an arbitrary initial configuration, the system reaches af- 
ter a relatively short transient time a state of the type 
. . . AABBCCAAAB ... in which A, B and C domains 
are located to the right of C, A and B domains, respec- 
tively. Due to the bias q, the domain walls . . . AB . . ., 
. . . BC . . ., and . . . CA . . ., are stable, and configurations 
of this type are long lived. In fact, the domains in these 
configurations diffuse into each other and coarsen on a 
time scale of the order of g^^, where L is a typical do- 
main size in the system. This leads to the growth of 
the typical domain size as (lni)/|lng|. Eventually the 
system phase separates into three domains of the differ- 
ent species of the form A . . . AB . . . BC . . .C. A finite 
system does not stay in such a state indefinitely. For 
example, the A domain breaks up into smaller domains 
in a time of order q-^i-n{NB,Nc} _ ^j^g thermodynamic 
limit, however, when the density of each type of parti- 
cle is non vanishing, the time scale for the break up of 
extensive domains diverges and we expect the system to 
phase separate. Generically the system supports particle 
currents in the steady state. This can be seen by con- 
sidering, say, the A domain in the phase separated state. 
The rates at which an A particle traverses a B (C) do- 
main to the right (left) is of the order of q^^ (9^^)- The 
net current is then of the order of q^'^ — q^'^ , vanishing 
exponentially with N . This simple argument suggests 
that for the special case Na = Nb = Nc the current is 
zero for any system size. 

In the following, we show that the dynamics of the 
model satisfies detailed balance for the special case Na — 
Nb — Nc- The Hamiltonian governing the steady state 
distribution is found to have long range interactions. It 
is demonstrated that phase separation takes place in the 
thermodynamic limit. The simple considerations given 
above indicate that phase separation exists even when 
the number of particles of the three species are not equal. 

B. Special case, Na = Nb — Nc : We specify a config- 
uration by the set of numbers {Xi} — {Ai, Bi,Ci}, where 
Ai , Bi and Ci are equal to one if site i is occupied by par- 
ticle A, B or C respectively and zero otherwise. We show 
that the dynamics satisfies detailed balance with respect 
to the Hamiltonian Ti. given by: 

JV-l N 

ni{X,}) =Y. [C*'^^ - C-.A, + B,A,]. (2) 

i=l j=i+l 

So that, in the steady state the probability of finding the 
system in a configuration {Xi} is given by: 

M/^.({X,}) = Z^lg^«^•», (3) 

where the partition sum is given by Zjv = 
"^{x } g'^^''-'^'^-' and the sum runs over all states in which 
A^^ = Nb — Nc- Although expression (||) for the Hamil- 
tonian is not manifestly translationally invariant, a form 
that is clearly invariant under translations can be simply 
obtained from (0) but is less convenient for our purposes. 



Before proving Eqs.(||,||) let us make a few comments. 
The A^-fold degenerate zero energy ground states of the 
Hamiltonian (|^) comprise the fully separated configu- 
ration A - . - AB - - - BC - - .C and any translation of this 
configuration. Starting from a ground state, any other 
configuration may be obtained by successive permuta- 
tions of nearest neighbor particles. The energy of such 
configuration may be calculated by noting that any par- 
ticle exchange against the bias costs one unit of energy 
while an exchange in the direction of the bias results in 
a gain of one unit. The maximal energy of this Hamilto- 
nian is A^^/9 and corresponds to a fully separated state 
A. - - AC . - - CB - - - B. It is of interest to note that al- 
though the dynamics defined in (Q) is local, the resulting 
Hamiltonian (H) is long range, in which each particle in- 
teracts with all other particles. Also the interactions are 
asymmetric in the sense that the Hamiltonian is not in- 
variant under space inversion. The temperature of the 
system is given by T — 1/| lng|. 

To prove Eqs.(||,^) one can check that WN{{Xi}) sat- 
isfies detailed balance with respect to the dynamics (|^). 
For example it is easy to see that an exchange of par- 
ticles in the bulk of the lattice . . . AB ... ... BA . . . 

satisfies qWNi---AB.-.) — Wn{---BA.--)- Similarly, parti- 
cle exchange between site 1 and N, say, A - - - C ^ C - - - A 

satisfies qWNiA C) = WnIC A), as long as Na — 

Nb — Nc- Note that when the three densities are not 
equal, the system supports a current, and thus it cannot 
satisfy detailed balance. 

In order to demonstrate phase separation it is suffi- 
cient to show that for large finite k the two point density 
correlation function satisfies 

\im {{A,Ak)~{Ai){Ak))>0- (4) 

N—*-oo 

In our case this implies limAr^oo(^i^fe) > 1/9 for large 
k- In fact we show below that in the limit of large N and 
for any finite k: 

{A,Ak) ^l/3-0{l/N), (5) 

demonstrating the existence of phase separation in the 
model. This result indicates that the phase separation is 
complete, in the sense that the probability of finding, say, 
a _B or a C particle a large distance inside the A domain 
vanishes in the thermodynamic limit. 

We proceed by first calculating the partition sum, Zn, 
showing that for q < I and to leading order in A^ 

Zm - N/ [iq)^f , (6a) 

where 

{q)oo = lim (l-g)(l-g2)...(l_g»). (6b) 

n — >oci 

The partition sum for g > 1 is obtained by replacing q 
by 1/q in this expression. Note that the partition sum 
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FIG. 1. The I — 5 ground state for an A'' = 21 system. 

is proportional to N rather than being exponential in 
N. This is a consequence of the long range interactions; 
low energy excitations are localized about the domain 
boundaries and their degeneracy does not increase with 
system size (in contrast to short range models). At q — 1, 
however, all configurations are equally probable and the 
partition sum is exponential in N. 

The calculation of the partition sum is greatly sim- 
plified by noting that configurations with energy larger 
than aN, where a > is a constant, can be neglected 
in the thermodynamic limit. Here we show that this is 
the case for q < (l/3)l/^ A more detailed analysis com- 
pletes the proof for all q ^ I |jl^. Consider the sum 

Zm>aN = Sm=aAr-(-i -^l™)?"' where D{m) is the num- 
ber of states of energy to. Noting that D{m) is bounded 
from above by the total number states which in turn is 
bounded by 3^, one finds that Zm>aN decays exponen- 
tially with N as long as q < (1/3)^/". Hence, the parti- 
tion sum may be replaced by the truncated one as long 
as exponentially small corrections in N are of no interest. 

To calculate the truncated partition sum wc note that 
states with energy m < N/3~ 1, may be uniquely decom- 
posed into TV disjoint sets, where each set corresponds to 
one of the N ground states of the system. Any state 
within a given set may be obtained from the correspond- 
ing ground state by successive permutations of nearest 
neighbor particles, in a way that the energy always in- 
creases along the path of intermediate states. 

The sets are labeled hy I = 1 . . . N, the position of the 
rightmost A particle in the A domain in the ground state 
configuration of that set (Fig.(l)). Due to the transla- 
tional invariance of the model the partition sum can be 
written as Zn = NZ -|- e~^'^^^ where Z is the parti- 
tion sum obtained by summing over one of the N sets of 
configurations. We now turn to calculating Z. Consider 
an excitation of energy m created only at one domain 
boundary, say AB. The excitation can be formed by one 
or more B particles moving into the A domain (equiva- 
lently A particles moving into the B domain). A moving 
B particle may be considered as a walker. The energy 
of the system increases linearly with the distance trav- 
eled by the walker inside the A domain. An excitation 
of energy m at the AB boundary is formed by j walkers 
passing a total distance of to. Hence, the total num- 
ber of states of energy to at the AB boundary is equal 
to the number of ways P{m) of partitioning an integer 
TO into a sum of (positive) integers. When considering 
all three interfaces in the ground state, the number of 



states composed of excitations of energy to^, z = 1, 2, 3 at 
the three interfaces is given by P(mi)P[m2)P{rnz)- Note 
that since only excitations with total energy smaller than 
A''/3 — 1 are considered, a walker can not travel from one 
boundary to another, and local excitations at the three 
boundaries are independent. The partition sum Z, after 
taking the thermodynamic limit, is given by 

oo m 

Z=Y. 9" H P(toi)P(to2)P(to3)5 (7) 

m— rrii—Q 

This sum may be rewritten as 

oo 

Z^{Y^q^P{m)f. (8) 

m=0 

Using the classical result for the generating function of 
P{m) 1^, attributed to Euler, we obtain Eqs.(||). Note 
that excitations of energy up to N/S — 2 are properly ac- 
counted for. Once excitations of greater energy are con- 
sidered, triplets of the form . . . ABC . . . can be formed. 
These triplets can be moved around the lattice with no 
energy cost so that the degeneracy is not accounted for 
correctly in Eqs.(|^). However, our results above imply 
that these contributions are unimportant for q < (1/3)"^. 
A more detailed analysis |l^ confirms that for all q < 1, 
Eqs. (|^) give the partition sum up to corrections expo- 
nentially small in N. 

In order to calculate the correlation function consid- 
ered in the next paragraphs, it is useful to introduce a 
partial sum Zs. It is defined as the partition sum under 
the constraint that one of the walkers at the AB interface 
has traveled a distance s. One can show that 

Zs = q'Z. (9) 

We now turn to the correlation function (AiAk) = 1/3 — 
{AiBk) — (AiCk)- In order to prove Eq. (^) we show 
that the correlation functions {AiBk) and {AiCk) are 
of 0{\/N). This is done explicitly for {AiBk). Similar 
considerations yield the same result for {AiCk)- Define 
{AiBkji as the correlation function calculated within the 
set of states related to the ground state I. Neglecting 
corrections exponentially small in N we can write 

1 ^ 

{AiB,) = -Y^{A,B,)i. (10) 
1=1 

We now outline the proof that X^z^i (^i^*:)' finite 
so that for large N , (AiBk) vanishes like 1/A^. Consider 
first the 1 = 1, Af/3 terms in the sum ( p^ ) for which 
Ai = 1 in the ground state. The sum of these terms may 
be broken into two parts in the following way: 

N/3 k-1 N/3 

Y,{AiBk)i = Y,iAiBk)i + Y,iAiBk)i. (11) 

1=1 1=1 l=k 
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To bound this sum we note that the first term con- 
tains contributions from states corresponding to ground 
states in which Ai ^ Bk — \. Here we use the bound 
{AiBkji < 1 for I = 1, ...k — 1. The second term contains 
contributions from states corresponding to ground states 
in which Ai = 1 but = 0. Here, only excited states 
in which one of the B walkers travels at least a distance 
I — k + 1 into the A domain may contribute to {AiBk)i. 
Thus, for A; < ^ < N/3 the correlation function {AiBk)i 
may be bounded from above by 'YlT=i-k+i -^s/ Z. Com- 
bining the bounds for the two parts one finds 



Y,{AiBk)i <k 

1=1 



N/3 oc 

E E 

l=k s=l-k+l 



(12) 



The geometrical series in (|l^ ) converges for q < 1 and 
thus the sum '^f^i {■^iBk)i is bounded in the thermo- 
dynamic limit. Similar considerations yield a bound 
for the remaining 2N/3 terms in ( |lO| ) demonstrating 
that {AiBk) 0(1/7V). We have thus shown that 
(AiAk) = 1/3 - 0(1/A^) for large TV and any finite fc, 
demonstrating phase separation in the model. 

So far we have presented some results and bounds in 
the thermodynamic limit. It is also of interest to ob- 
tain exact results for finite rings to investigate how the 
limit is approached. We do this by applying a matrix 
ansatz method which has recently been introduced for 
studying Id non-equilibrium systems |11|. Generaliz- 
ing this approach to replace the matrix product used 
as steady state ansatz by a tensor product we found 
that we could apply the method to the model for 
the case Na = Nb = Nc- Details of these calculations 
will be presented elsewhere jl^ , however we would like to 
present some exact results for small systems obtained by 
this method. We have calculated the correlation function 
{AiApf/2) which provides a measure of phase separation. 
In a disordered state this correlation function is equal to 
(iV — 3)/9(iV— 1) and approaches 1/9 in the large N limit. 
It should be smaller for a phase separated state. In fact 
we find it approaches zero over a range of q values which 
increases with N (Fig. (2), inset). To investigate the fi- 
nite size scaling near the q — 1 (infinite temperature) 
critical point it seems natural to choose Nhiq as a scal- 
ing variable. This variable represents the ratio of domain 
wall width (1/| lnq|) to domain size {N/3). In Fig. (2) the 
scaling collapse for small systems is illustrated. 

C. Discussion: The analysis presented above dealt 
with the case Na = Nb = Nc- In the general case where 
the densities of the three species are not equal, detailed 
balance is not satisfied. However, the heuristic arguments 
for phase separation given at the beginning of the Let- 
ter are expected to hold for the general case provided 
that none of the densities is zero. Namely, configurations 
of the type . . . AAABBBBCC'AAA ... are stable, and 
the time for a totally phase separated state to break up 
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FIG. 2. The correlation function C = (yliAjv/2) as a 
function of the scaled variable N\a.q for A'^ = 30, 36, 42, . . . , 84. 
The inset shows the same data plotted against q. Note that 
at g = 1 the curves approach 1/9 as A'' increases (see text). 



grows exponential in N . Numerical simulations of the 
model supports the existence of phase separation | |To|] . 
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